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Abstract. We investigate the rapidity and transverse momentum distributions of 
the prompt photon production at the CERN LHC energies considering the current 
perturbative QCD approaches for this scattering process. Namely, we compare 
the predictions from the usual NLO pQCD calculations to the the color dipole 
formalism, using distinct dipole cross sections. Special attention is paid to parton 
saturation models at high energies, which are expected to be important at the 
forward rapidities in pp collisions (y^ — 14 TeV) at the LHC. 

1 Introduction 

There are recent experiments collecting prompt photon data, as the PHENIX collaboration at 
^/s = 200 GeV at RHIC (collecting data both for the inclusive [1] and the isolated case |2]) and 
the D0 collaboration [3] {^/s = 1.96 GeV) in the Tevatron (measuring isolated prompt photons 
whose transverse momenta px range from 23 to about 300 GeV) . A photon is said to be prompt 
if it does not originate from the decay of a hadron, such as a tt" or rj, itself produced with large 
transverse momentum. From theoretical point of view, the cross sections for producing such 
isolated photons have been proven to still fulfill the factorisation property, and are finite to all 
orders in perturbation theory. The study of prompt photons carrying large values of transverse 
momentum has long history and it is related to deep-inelastic-lepton scattering (DIS), Drell- 
Yan process production and jet production as an important probe of short-distance hadron 
dynamics. The great appeal of prompt photons is that they are point-like, colorless probes of 
the dynamics of quarks and giuons, once that escape unscathed through the colored medium 
of the high-energy collision. Thus, they can be a powerful probe of the initial state of matter 
created in heavy ion collisions, since they interact with the medium only electromagnetically 
and therefore provide a baseline for the interpretation of jet-quenching models. Prompt photons 
offer access to the spin-dependent and spin-averaged gluon densities of hadrons, where one of 
the two leading-order partonic direct subprocesses feeds directly from the gluon parton density 
through the 'Compton' subprocess, qg jq. It should be stressed that the simplicity of prompt 
photons is compromised with the fragmentation contributions, where fragmentation is a long- 
distance process in which a hard photon brems off a final-state quark (or gluon). The photon 
emerges as part of a jet if the opening angle between the quark and photon is too small and 
this 'showering' process is parameterized by a non-perturbative single-photon fragmentation 
function D^{z, /if)- 
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Calculations exist of prompt single photon production at next-to-leading order (NLO) in 
QCD for both the direct and the fragmentation contributions, as the approach is encoded in 
the Monte Carlo program JETPHOX [4,5,6 . NLO pQCD calculation provides a fairly good 
description of the world-data on px > 5 GeV hadroproduction spectra [4], where the mid- 
rapidity cross section is sensitive to parton distributions at x ~ 2pT/y/s and ~ p'^ and 
at a given Q^, other values of x can be probed by measuring photons at non-zero rapidity. 
Accordingly, the gluon densities in a proton and in a nucleus are fundamental ingredients in 
order to compute the corresponding hard-process observables in pp, pA and A A collisions. The 
gluon distribution xg(x, Q^) is fairly well known within a few percent accuracy in the range 
X ^ 10~^-10~^ and ~ 10-10^ GeV^ T, which is precisely the kinematical domain covered 
by most hard processes at the LHC. However, very little is known on the nuclear gluon density 
per nucleon, xgAix, Q^) [8]. One way to better constrain the nuclear gluon pdf is though prompt 
photon production. In Ref. [9] the prompt photon production in p-Ph collisions at LHC and 
d-Au collisions at RHIC has been explored, using that the nuclear production ratio of isolated 
photons in p-A collisions can be simply approximated as a linear combination of gluon pdf 's 
and structure functions in the nucleus A over those in a proton. Performing the calculation in 
pQCD at NLO they have checked that such an approximation is correct up to a few-percent 
accuracy and makes this observable an ideal tool to measure gluon shadowing at RHIC and 
at LHC. A LO QCD study of prompt photon production in pA and AA collisions including 
more recent nuclear pdf's is done in ref. [10] . and it is shown that the amount of shadowing 
and antishadowing is very distinct for the various nuclear pdf's available, so these processes are 
useful! to discriminate among the different nuclear pdf's. Moreover, a precise knowledge of the 
gluon pdf's is crucial for a better understanding of the QCD evolution, mostly the possible non- 
linear QCD evolution at small x and small , where the gluon density in the nucleus becomes 
large and starts to saturate at momentum scales near Qsat (the so-called "saturation scale" ) [11] . 
The influence of nonlinear gluon evolution in prompt photon production is estimated in Ref. 
|12j . where we find an enhancement of the prompt photon cross section for LHC energies. 

Recently, it has been verified that the color dipole approach can successfully describe inclu- 
sive photon production in pp collisions at midrapidities. In Ref. "13] it was showed that both 
direct photon production and Drell-Yan dilepton pair production processes can be described 
within the same color dipole approach without any free parameters. Such a formalism, de- 
veloped in [14] for the case of the total and diffractive cross sections, can be also applied to 
radiation |15j . Although in the process of electromagnetic bremsstrahlung by a quark no real 
quark dipole participates, the cross section can be expressed via the more elementary cross 
section adip of interaction of a QQ dipole. It was also showed that the color dipole formulation 
coupled to the DGLAP evolution provides a better description of data at large transverse mo- 
mentum compared to simple saturation dipole models. In contrast to the parton model, any 
photon fragmentation function was included for computing the cross section, since the dipole 
formulation already incorporates all perturbative (via Pomeron exchange) and non-perturbative 
radiation contributions. Using similar approach, i.e. the Color Glass Condensate (CGC) for- 
malism, in Ref. |16j the authors presented the predictions for prompt photon production using 
a model for the scattering dipole amplitude which describes quite well the hadron production. 
As stressed in Refs. [T7 ] fT8 ] fT9 ] . electromagnetic probes of the CGC are crucial to determine the 
dominant physics in the forward region at RHIC and LHC. They estimated the ratio R^a for 
photon production at forward rapidities for RHIC and LHC energies and compare its behavior 
with that predicted for hadrons and predicted the photon to pion production ratio and study 
its pt dependence. Therefore, such an approach is quite promissing to predict observables to 
measured at the LHC. 

Our goal in this contribution is twofold. Firstly, we perform a comparison between the 
NLO QCD approach to the color dipole formalism pointing out the possible deviations and 
their origin. Secondly, we investigate in detail the color dipole approach using distinct imple- 
mentations for the color dipole cross section, discussing several phenomenological aspects. The 
paper is organized as follows. In next section we summarize the main formula for both NLO 
pQCD approach and color dipole formalism. In last section we show our numerical results and 
predictions. Finally, we summarize our main conclusions. 
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2 Prompt photon production 

Let us summarize the theoretical treatment for prompt photon production in pp collisions in 
high energy colliders. We start with the usual parton model description of underlying processes 
and after that we present the alternative color dipole formalism. We call attention that the 
dipole approach is well suited for high-energy processes, i.e. small X2, and its range of validity 
is expected to be near X2 < 0.1. The main advantage in such an approach is to describe 
simultaneously the direct photon and dilepton production in the same framework. 



2.1 NLO QCD approach 



Schematically, the production of a prompt photon proceeds via two mechanisms. The first one 
is the 'direct' mechanism, where the photon behaves as a high px colourless parton, i.e. it takes 
part in the hard subprocess, and it is most likely to be well separated from any hadronic envi- 
ronment. The second one is called 'fragmentation' mechanism, where now the photon behaves 
as a kind of (anomalous) hadron, i.e. it results from the coUinear fragmentation of a coloured 
high pt parton, and is it most probably accompanied by hadrons - unless the photon carries 
away most of the transverse momentum of the fragmenting parton. The fragmentation contri- 
bution emerges from the calculation of the higher order corrections to direct mechanism in the 
perturbative expansion in powers of the strong coupling ag. At higher orders, final state multi- 
ple coUinear singularities appear in any subprocess where a high pr parton of species k (quark 
or gluon) undergoes a cascade of successive coUinear splittings ending up with a splitting into 
a photon. These singularities are factorised to all orders in as according to the factorisation 
theorem, and absorbed into fragmentation functions of parton fc to a photon, Dry/^iz, M^), 
defined in some arbitrary fragmentation scheme, at some arbitrary fragmentation scale . 
The point-like coupling of the photon to quarks is responsible for the well-known anomalous 
behaviour of D^/i,{z, M^), roughly as aem-ce7^{Mp) when the fragmentation scale M^, chosen 
of the order of a hard scale of the subprocess, is large compared to 0(1 GeV). The differential 
cross section in transverse momentum p^ and rapidity rj can thus be written as: 



a{pp -> 7 + X) = ifiR, M, 



,(frag) 



k=q,q,g 



{^iR,M,M^)®D^^kiM,) (1) 



where a 



(frag) 



describes the production of a parton k in a hard collision. The arbitrary pa- 



rameters fiR, M and are respectively the renormalisation, initial-state factorisation, and 
fragmentation scales. The dependence of the NLO predictions with respect to hr, M and 
is simplified if we take all scales to be equal and they will be noted ii. We rely on the calculation 
of both direct and fragmentation at next-to- leading order (NLO) accuracy POIIS]. which takes 
the form (77 is the photon rapidity) 



da da(dir) ^^(frag) 



dpj.dri 



dpj.dri dpj,dr] 



(2) 



where 



dp^di] 



E 

daij 
dp^drj 



dxidx2 Fi/hJxi,M) Fj/h^{x2,M) 



2tt 



2tt 



(3) 



and 



da(frag) 
dp^drj 



E 

l,j,k=q,q,g 

ddh 



dz 



dxidx2— F,/h,{xi,M) Fj/h^{x2,M) D^/kiz,M^) 



27r 



dp^drj 



2tt 



(4) 
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where Fi/^i 2 (2;, M) are the parton distribution functions of parton species i inside the incoming 
hadrons /ii,2, at momentum fraction x and factorisation scale M; asinn) is the strong coupling 
defined in the MS renormalisation scheme at the renomalisation scale fin. The knowledge of 
A , e.g. from deep-inelastic scattering experiments, completely specifies the NLO expres- 

M3 

sion of the running coupling as{fiR). The NLO correction terms to direct and fragmentation, 

( dir ) ( ircLff ) 

Kl^ ' [U and Kl^P' M respectively, are known and their expressions in the MS scheme 

will be used. In the calculation of numerical results in next section, we will show the direct 

and fragmentation contribution, presenting the uncertainties associated with different choices 

for the scale fi. 



2.2 Color dipole formalism 

The color dipole formalism can be used to compute direct photon production as it is also 
applied to describe radiation processes [15] . The transverse momentum px distribution of pho- 
ton bremsstrahlung in quark- nucleon interactions, integrated over the final quark transverse 
momentum, was derived in 23] in terms of the dipole formalism. 



dilna)d^PT -{2n) ^^^^^ 



X ^{'^dip{x, ari) + adip{x, ara)} - ^(Tdip{x, a{ri - ra)), (5) 

where ri and r2 are the quark-photon transverse separations in the two radiation amplitudes 
contributing to the cross section, (Tdip. The parameter a is the relative fraction of the quark 
momentum carried by the photon, and is the same in both amplitudes, since the interaction 
does not change the sharing of longitudinal momentum. In the equation above, T stands for 
transverse and L for longitudinal photons. The energy dependence of the dipole cross section, 
which comes via the variable x — 2(pi ■ q)/s = (jpT / \fs) e.~y , where pi is the projectile four- 
momentum and q is the four-momentum of the dilepton, is generated by additional radiation 
of gluons which can be resummed in the leading ln(l/2:) approximation. 

In Eq. ([5]) the light-cone wavefunction of the projectile quark 7(7 fluctuation has been de- 
composed into transverse (/)^g(a, r) and longitudinal (t>^^{a, r) components, and an average over 
the initial quark polarization and sum over all final polarization states of quark and photon is 
performed. For direct photons, only transverse component contributes and this wavefunction 
component (j)'^^{a, r) can be represented at the lowest order as: 

J2 <,^(a,ri)0^,(a,r2) = |^ lmyKo{en)Ko{er2) + [1 + (1 - af]e^'^K^{en)K,ier2)] 

inj 

in terms of transverse separation r between photon 7 and quark q and the relative fraction a of 
the quark momentum carried by the photon. Here Ko,i{x) denotes the modified Bessel function 
of the second kind and it has been introduced the auxiliary variable = a^m^, where ruq is 
an effective quark mass. Here, we take = 0.2 GeV for direct photon production. 

The hadron cross section can be obtained from the elementary partonic cross section Eq. ^ 
summing up the contributions from quarks and antiquarks weighted with the corresponding 
parton distribution functions (PDFs), 

dajpp jX) ^ '^pp^p_ ^2^ d(T'i^{q qj) 
dyd'^PT J XX ct ^ a ' d{lna) d'^px ' 

where the PDFs of the projectile have entered in a combination which can be written in terms 
of proton structure function For the hard scale /z entering in the proton structure function 
in Eq. ([5]), we take = p"^ and the energy scale x of the dipole cross section entered in Eq. ^ 
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was set X = X2- For the proton structure function in Eqs. ([6]) we have taken the recent ALLM 
parametrization [24] , which is valid in the kinematic range we are interested in. The sensitivity 
to a different choice for F2 is very small. 

A crucial ingredient in the color dipole calculations is the dipole cross section. It is theo- 
retically unknown, although several parametrizations have been proposed. For our purposes, 
here we consider some analytical parametrizations. The simplest parametrizations rely on the 
geometric scaling property. That is, the dipole cross section is a function of a scaling variable 
TQsat{x) where Qsat is the so called saturation scale. It defines the transverse momentum scale 
where parton recombination physics is relevant and in general is modeled as Qsat oc x^^/^ (it 
grows with energy). A common feature on these models is that for decreasing x, the dipole cross 
section saturates for smaller dipole sizes, and that at small r, as perturbative QCD implies, 
(T ~ vanishes, i.e. the color transparency phenomenon 14 . In a general form, the dipole 
cross section can be parametrized as. 



(Tdip(x,Y\^) = (To 



1 — exp 



lL,{x) = {^)\ (7) 



where the quantity 7eff is the so-called effective anomalous dimension. For the celebrated GBW 
parametrization [25], 7eff = 1 and the remaining parameters are fitted to DIS HERA data 
at small x. This parametrization gives a quite good description of DIS data at x < 10~^. 
A very recent adjust using GBW model for 0.25 < Q'^ < 45 GeV^ and using quark 
masses m, = 0.14 GeV and rric = 1.4 GeV gives ctq = 23.9 mb, = 1.11 x 10"^, A = 0.287 
(xVdof = 1.58). 

The main difference among the distinct phenomenological models using parametrizations as 
Eq. ([7]) comes from the predicted behavior for the anomalous dimension, which determines the 
transition from the non-linear to the extended geometric scaling regimes, as well as from the 
extended geometric scaling to the DGLAP regime. It is the behavior of 7 that determines the 
fall off with increasing pT of the hadronic cross section. The current models in the literature 
consider the general form 7efi = 7sot + ^(2^, r; p^) , where ^sat is the anomalous dimension at the 
saturation scale and A mimics the onset of the geometric scaling region and DGLAP regime. 
One of the basic differences between these models is associated to the behavior predicted for A. 
While the models proposed in Refs. [27,28,29] assume that A depends on terms which violate 
the geometric scaling, i.e. depends separately on r and rapidity Y — ln{l/x), the model recently 
proposed in Ref. [3^ by D. Boer et al. consider that it is a function of rQsat- 

7eff = Jsat + (1 - Jsat) _ X) + fo ' 

where u) = pr/Qsat and the two free parameters a = 2.82 and h ~ 168 are fitted in order 
do describe the RHIC data on hadron production. The remaining parameters are taken from 
the GBW parametrization. Thus, for large pr the Boer at al. dipole cross section reproduces 
the GBW model. They differ at small pT < Qsat compared to the saturation scale, where 

7eff ~ Isat ^ 0.63. 

In order to investigate the effect of QCD evolution in the dipole cross section we also will 
include in our studies the impact parameter saturation model 26]. It has been successful in 
describing the exclusive vector meson production. The dipole cross section in this model is given 
by: 



crdip{x,r) =2 j '^^^ 



^2 

1 - exp ( -—r'^as{p^)xg{x, n^)T{h) 



(9) 



where the scale fj? is related to the dipole size r by — 4/r'^ + fiQ. The gluon density, xg(x, jj.'^), 
is evolved from a scale /Iq up to fj,^ using LO DGLAP evolution without quarks. The initial 
gluon density at the scale /ip is taken in the form xg{x, /ip) — Ag x~^o (1 — x)^'^ . The values of 
the parameters /Zq, Ag, and Xg are determined from a fit to F2 data. For the light quarks, the 
gluon density is evaluated aX x = xb (Bjorken-x), while for charm quarks, a; = Xb^^ + ^to^/ Q^). 
The LO formula for the running strong coupling as{lJ^) is used, with three fixed flavours and 
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PTlGeV] IGeVJ 
(a) (b) 

Fig. 1. (a) The invariant cross section for prompt photon production in pp collisions at RHIC compared 
to NLO QCD and color dipole (using GBW model) results, (b) Predictions for the differential cross 
section d^a{pp — > j X)/dpT d-q using the NLO QCD and color dipole formalism at the Tevatron (|7;| < 
0.9). 



^QCD — 0.2 GeV. The proton shape function T{b) is normalised so that J (PhT{h) = 1 and 

a Gaussian form for T(6) is considered, Tc{h) ~ ^^^^ c . A recent adjust using b-SAT 
pS] model for 0.25 < < 650 GeV^ using quark masses niq = 0.14 GeV and rric = 1-4 
GeV gives = 1-17 GeV^, Ag — 2.55 and Xg = 0.020, producing a good quality ajust with 
X^/dof = 1.21. In next section we will use the phenomenological models presented here to 
compute the px spectrum of direct photon production in pp collisions. 



3 Numerical results 

Let us now present some numerical results concerning direct photon production in central 
rapidities. In Fig. [T]-a, we show the result for the direct photon production in pp collisions 
at RHIC obtained by the NLO QCD approach and color dipole picture. The invariant cross 
section, Ed^a/(PP, is computed as a function of transverse momentum for central rapidity, 
y = 0, at ^/s = 200 GeV and data are taken from Ref. [3T]. For the NLO QCD calculation, we 
test the theoretical uncertainty using two distinc scales: fi = pt/2 (dotted line) and = px 
(solid line). We have included also the photon fragmentation contribution [32j with scale /i = px 
(long dashed line). In these calculations, we used the CTEQ5M parton distribution functions 
(small changes using alternative PDFs sets). We verify relatively small uncertainties coming 
from set of PDFs and/or hard scale. For the color dipole calculation, we reproduce the results 
on Ref. p3] where the GBW dipole cross section is considered (dot-dashed hue) . The deviation 
between NLO QCD and color dipole starts to be sizable at transverse momenta of order 10 
GeV, reaching one order of magnitude at pr = 25 GeV. In Fig.[T]-b, the differential cross section, 
(Pa/drjdpT, is computed for pp collisions at the Tevatron, y/s = 1.8 TeV and \ri\ < 0.9, for 
different choices of scale and parton distributions. The notation is the same as in the previous 
figure. The experimental points are CDF data [33^. The NLO QCD describes the data quite 
well, while the color dipole (GBW without evolution and with threshold correction) result seems 
to give an incorrect behavior. Next we show how this problem can be cured within the color 
dipole picture. 

In Fig. [2]-a we show the invariant cross section (central pseudorapidity 77 = at RHIC 
[31] , ^/s = 200 GeV) obtained from several implementations of the dipole cross section taken 
from recent phenomenological works. For instance, we the following models: (i) the baseline 
calculation of Ref. [13] (Kopehovich et al., dotted solid line); (ii) the new fit with GBW dipole 
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Pt [GeV] p, [GeV] 

(a) (b) 

Fig. 2. (a) The invariant cross section EdPaijjp 'yX)/d?P for pp collisions at RHIC compared 
to different implementations of dipole cross sections within the color dipole picture (see text), (b) 
Predictions for the differential cross section cfa{pp j X)/dr]dpT using the several implementations 
for the dipole cross section for Tevatron energy. 



cross section, corrected by large-x threshold (dot-dashed line); (iii) the impact parameter sat- 
uration model which includes QCD evolution in the dipole cross section, labeled here b-SAT 
(solid line); (iv) the D. Boer et al. model (labeled BOER et al., dashed line) which uses a 
running anomalous dimension. The main deviation among the models occurs at large pt- We 
first verify that the supposed incorrect behavior of GBW model discussed in Ref. [I3] is due 
to the lack of a large- a; threshold corrections. We have multiplied the GBW dipole cross sec- 
tion by a factor (1 — ^2)^ and this fact produces large suprression of cross section at large 
Pt- For sake of ilustration, we show the vertical line at pr ~ 2 GeV which corresponds to 
X2 — Pt/\^ = 10^^, that is the expected limit of the color dipole approach. The RHIC data 
are all clearly within the large-x region and, therefore, a threshold correction should be added in 
a consistent phenomenology. The BOER et al. model improves the data description al small px 
but is similar to GBW towards large transverse momenta. The reason for that is that at large 
Pt the effective anomalous dimension in both models coincide, giving 7=1. The b-SAT model 
gives similar results as GBW at small px, wheras the data description is improved at large px- 
This property is already know in studies from Ref. |13| . where the QCD evolution present in 
BGBK saturation model gave better data description. Finally, we show that saturation physics 
is not directly relevant for RHIC at midrapidity by writing down the values of the saturation 
scale in the kinematic range of data. We consider for simplicity the saturation scale from GBW 
mode, Qlat = {xo/x2)^ = {xq^/s /pr)^ . We got 0.1 < Q^^t < 0.5 GeV^, which is very smaU 
compared to the transverse momenta 4 < < 100 GeV^, Therefore, saturation effects do not 
play an important role at RHIC midrapidity. 

In Fig. [H-b, we performed the same study for Tevatron. It is shown the differential cross 
section cfa/dr]dpT as a function of transverse momentum for CDF energies ^/s = 1.8 TeV |33) . 
At lower transverse momentum pT < 40 GeV the GBW dipole model can reproduce rather 
fairly the experimental data, and at higher px values DGLAP present in the b-SAT evolution 
significantly improves the results. It should be stressed that even at Tevatron energies, the large- 
X effects are still important. We show the vertical line at px — 20 GeV which corresponds to 
X2 = 10^^, telling us that the large px data at Tevatron are out of small- a; region. Accordingly, 
we have used the dipole models corrected by threshold factor as referred before. The b-SAT 
model already contains such a correction by construction. The same statements about the role 
played by saturation efects remains valid for Tevatron at midrapidty, where the saturation 
scale is in the range 0.2 < Ql^t ^ 0-8 GeV^ . As a partial conclusion, at low transverse 
momentum px < 10 GeV saturation models predictions are almost identical, but at higher px 
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the dipole parametrization including DGLAP evolution bends down towards the experimental 
points improving the result. 

In Fig.[3]-a, the differential cross section, (Pa/drjdpT, is computed for collisions at central 
rapidity at the LHC, ^/s = 14 TeV, for the lower pT region up to 30 GeV. The notation is the 
same as in figure [TJ It is verified that the sensivity on the scale for the NLO QCD calculation 
is larger than for RHIC and Tevatron, mostly at smaller pT- The deviation between QCD and 
color dipole remains the same at the LHC. In Fig. [3]-b we show the dipole models predictions 
(GBW and b-SAT) for inclusive prompt-photon production at midrapidities, and for LHC 
energies y/s — 14 TeV. At lower transverse momentum px < 30 GeV the GBW dipole model 
produces similar results as the b-SAT model and at higher px values DGLAP evolution becomes 
significantly important. It is verified that even at midrapidities at the LHC, the X2 values are 
small as it can be shown by the vertical line at px — 130 GeV. Once again, despite the saturation 
scale reaching values above a few GeV at central rapidities at the LHC, saturation effects should 
be negligible. The correct place to search such effects is looking at very forward rapitidies in 
the pp mode or in the nuclear mode, where saturation scale should be enhanced by a factor 
A^^"^ (a factor 6 or 7 for Lead). 

Finally, we will investigate semi-analytical calculations which are allowed in color dipole 
picture in the color transparency region. First, we can write the explicit expression for px 
distribution using Eq. ([6|) and the expressions for the transverse momentum px distribution of 
photon bremsstrahlung in quark- nucleon interactions, Eq. (O, and transverse light cone wave 
function, Eq. ([6]). It reads as: 



da (pp jX) Uem da 



dyd'^px 2 



2 4 



+ e'^) 4e 



+ [l + {l-af 



{Px + e') 



Ii ^ el2 



(10) 



where we recall the auxiliary function s — amq. The quantities Ti,2,3 are Hankel's integral 
transforms of order (I1.2) and order 1 (X3) given by: 



dr rJo (px r)Ko (e r) a^ip , ar) , 



drr^Jo{pxr)Ki{er) a^ip {x2 , ar ) , 



X3 = / drrJi{pxr)Ki{sr)adip{x2,ar). 
Jo 



(11) 
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If we consider that the current experimental resuhs for midrapidities select a kinematic 
interval where the dipole cross section is dominated by the color transparency region, we can 
use the GBW parametrization and take its small-r limit to compute analytically the integrals 
in Eq. (jlip . In this case, we can take the approximation a dip ~ o'o(»'^ Q^at) the region where 
Pt > Qsat and we get: 



^1 -'^Oy.at (p|_^£2)3' 



4e (£2 



Qlat 

VP' 



2pTe 



+ £2)3- 



(12) 



We have verified that analytical calculations of referred integrals can be also done if we 
consider a fixed 'effective' anomalous dimension. Taking the result presented in Eq. an 
explicit expression for the hadron cross section is obtained. 



da{pp jX) 
dycpTpT 



2 

sat 



27r2 



da 
a 



V a 



2 4 



Pt 



(p|+£2)4 



+ [1 + (1 - a)2 



Pt 



2(p2, + £2)4 



Therefore, the anomalous dimension for models relying on extended geometric scaling play 
an important role in the functional form of the pT spectrum of prompt photon production in 
hadron-hadron collisions. 

As a summary, we showed that direct photon production can be addressed in the NLO 
QCD approach and in the color dipole picture. In particular, direct photon production can 
be described within the color dipole approach without any free parameters by using dipole 
cross section determined from current phenomenology in DIS. In central rapidities at RHIC 
and Tevatron, saturation effects do not play a significant role for direct photon production 
at the given experimental range of pt- This situation can be changed at the LHC even at 
midrapidities. We also verified that the color dipole formulation including DGLAP evolution 
(b-SAT model) provides a better description of data at large transverse momentum compared 
to the dipole models using fixed or running anomalous dimension. We write down also analytical 
results relevant for the calculation of pt spectrum of direct photons. Our results corroborates 
the statements from Ref. |13J and call attention for the large-x effects that should be added to 
the corresponding phenomenology for RHIC and Tevatron. 
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